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ABSTRACT 

We report results of an extensive observational campaign of the 6 cm formalde- 
hyde maser in the young massive stellar object IRAS 18566+0408 (G37.55+0.20) 
conducted from 2002 to 2009. Using Arecibo, VLA, and GBT, we discovered 
quasi-periodic formaldehyde flares (P ~ 237 days). Based on Arecibo observa- 
tions, we also discovered correlated variability between formaldehyde (H2CO) and 
methanol (CH3OH) masers. The H2CO and CH3OH masers are not spatially co- 
incident, as demonstrated by different line velocities and high angular resolution 
MERLIN observations. The flares could be caused by variations in the infrared 
radiation field, possibly modulated by periodic accretion onto a young binary 
system. 
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Suhject headings: ISM: molecules - masers - radio lines - ISM: individual 
(IRAS 18566+0408) 

1. Introduction 

IRAS 18566+0408 (G37.55+0.20) is a massive star forming region located at a kinematic 
distance of 6.7 kpc (Araya et al. 2004). The region is visible in Spitzer GLIMPSE infrared 
data (Benjamin et al. 2003), and shows an outflow that has been detected in CO, SiO and 
in the 4.5 /im IRAC band (Araya et al. 2007c; Beuther et al. 2002; Zhang et al. 2007). 
At smaller scales. Very Large Array (VLA, NRAO)0 observations reveal a radio continuum 
source, possibly an ionized jet (Araya et al. 2007c). The bolometric luminosity of the object 
is ~ 6 X IO^'Lq (equivalent to an 08 ZAMS star; Zhang et al. 2007; Sridharan et al. 2002). 

IRAS 18566+0408 harbors one of the few known 6 cm formaldehyde (H2CO) masers in 
the Galaxy (Araya et al. 2008). The H2CO maser in IRAS 18566+0408 showed a strong 
flare in 2002 (Araya et al. 2007b). Here we report results of a monitoring program of the 
H2CO maser intended to investigate the nature of the flare. 

2. Observations 

2.1. Arecibo Observations 

Using the 305 m Arecibo Telescop^, we monitored the 6 cm II2CO maser {JxaKc = 
lio - 111, i^o = 4829.6594MHz) in IRAS 18566+0408. The pointing position was R.A. = 
18^59™09.98^ decl. = 04°12'15.6"(J2000; Araya et al. 2005). The monitoring program had 
three intervals: May 2006 to April 2007, January 2008 to May 2008, and October 2008 to 
November 2009. A total of 48 runs were conducted (Figure 1; see Table 1, online version). 

We observed in position switching mode; 5 minutes on-source integration per scan with 
typically one or two scans per run. Observations prior to June 2008 were conducted with 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated 
by Cornell University under a cooperative agreement with the National Science Foundation. 
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the Interim auto correlator, dual linear polarization, 1.56 MHz bandwidth and 2048 channels 
(0.047km per channel). The data reduction was done using the CLASS programjfl The 
WAPP spectrometer was used since October 2008, with 3.125 MHz bandwidth and 2048 
channels (0.095 km s~^ channel width). The calibration of the WAPP data was done in IDL 
using Arecibo-IDL routines. The data were smoothed to a channel width of 0.19 km s~^. We 
observed B1857+129 in most of the runs to check the pointing and measure the telescope 
gain. The pointing error was less than 12" in general, and the gain varied between 5 and 9 
K/Jy in the 48 runs. The system temperature was typically 26 K. The Arecibo half power 
beam width (HPBW) is ~ 1' at 4860 MHz. We detected no significant variability of the 
linewidth or peak velocity of the maser. 

The 6.7 GHz CH3OH line (J^ = 5i - 60A+, = 6668.5192 MHz) was also observed 
with integration times on-source between 1 and 2 minutes. The same Interim and WAPP 
configurations described above were used for the CH3OH observations. The spectra were 
smoothed to a channel width of 0.07km s~^ (see Figure 2). The Arecibo HPBW is ~ 0.72' at 
6600 MHz. A detailed description of all CH3OH observations together with results from an 
ongoing monitoring program will be reported in the future. 



2.2. VLA Observations 

We used the VLA to observe the 6 cm H2CO maser in IRAS 18566 +0408 at seven epochs 
after 2006 (see Table 1, online version). The pointing positions were R.A. = 18*^59^10. 10*^, 
decl. = 04°12'12.0"(J2000) for the A and B array observations, and R.A. = 18^59^10.00% 
decl. = 04°12'25.0"(J2000) for the C array observations. The seven observations, together 
with two previous VLA observations from the literature (Araya et al. 2007b), are shown 
in Figure 1. We used a 2IF mode with bandwidth of 1.56 MHz (97kms~^) and 255 chan- 
nels (6.104kHz, 0.38 km s~^ channel width). We observed 3C286 and 3C48 as flux density 
calibrators, with assumed flux densities of 7.52 Jy and 5.47 Jy, respectively. We observed 
J1824+107 for amplitude and phase calibration. We measured flux densities between 0.74 
and 0.79 Jy for J1824+107 in the seven runs, i.e., a dispersion of ~ 6% which includes flux 
density calibration uncertainties and intrinsic variability of the quasar. All data reduction 
was done with the NRAO software package AIPS following standard spectral-line proce- 
dures. The H2CO maser positions measured with the VLA (A configuration) in 2003 and 
2007 agree within 50 mas. 



^CLASS is part of the GILDAS software package developed by IRAM. 
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2.3. GBT Observations 

The Green Bank Telescope (GBT) was used on five dates in 2008 to observe the H2CO 
maser (see Table 1, online version). The GBT H2CO flux density measurements from 2008, 
together with prior GBT observations from the literature (Araya et al. 2007b) are shown in 
Figure 1. 

The quasar J1851+005 was observed for pointing and focus adjustments. The pointing 
corrections were less than 8". The observations of April 07 and 15, 2008, were done in 
standard position switching mode with integration times on-source between 5 to 10 min. 
The observations in May and October 2008 were conducted in frequency switching mode, 
with integration time on-source between 10 and 30 min. We used the GBT spectrometer, 
with a bandwidth of 12.5 MHz (775kms~^), 8192 channels (0.094 km s~^ channel width) and 
dual circular polarization. The system temperature was ~ 20 K in all runs. After checking 
for RFI and consistency between the two circular polarizations, the data were averaged and 
smoothed to a channel width of 0.38 km s~^. All data reduction was done in IDL using the 
GBTIDL procedures!^ 



2.4. MERLIN Observations 

High angular resolution observations of the 6.7 GHz CH3OH masers in IRAS 18566+0408 
were conducted with Multi-Element Radio Linked Interferometer Network (MERLIN)! on 
April 03, 04, and 20, 2008, to determine the location of the CH3OH masers with respect to the 
H2CO maser. The pointing position was R.A. = 18^59™09.98^ decl. = 04°12'15.6"(J2000). 
The CH3OH masers were observed with a bandwidth of 0.5 MHz and 255 channels (1.96 kHz, 
0.088 km s~^ channel width). Each observing run lasted approximately 10 hours. The quasar 
1904+013 was observed as complex gain calibrator with a broad bandwidth (13 MHz) during 
approximately 2 minutes for every ^ 7 minutes observation of the CH3OH masers. We 
observed 3C84 in the narrow (0.5 MHz) and wide (13 MHz) setups to determine the solutions 
needed to transfer the phase calibration from the wide to the narrow band. 3C84 was also 
used as bandpass and flux density calibrator (a flux density of 15.7 Jy was assumed). 

Poor phase tracking affected the observations of April 03 and 04, 2008, thus the astrom- 
etry of the MERLIN observations discussed in this paper comes from the phase referenced 



''http://gbtidl.nrao.edu/ 

^MERLIN is a National Facility operated by the University of Manchester at Jodrell Bank Observatory 
on behalf of STFC. 
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observations of April 20, 2008. To increase signal to noise, we combined all MERLIN data 
and used the position of the brightest CH3OH maser from the 2008 April 20 observations as 
initial model for self-calibration. The synthesized beam of the self-calibrated data is 73 x 
42 mas, 22° position angle; the final rms is ~ 12mJyb~^. The self-calibrated data were used 
to create Figure 3. Including uncertainties in the phase reference source position, the tele- 
scope positions, and phase transfer errors, we estimate that the total astrometric uncertainty 
of our MERLIN data is 35 mas. 

3. Results 

We detected five new flare events and discovered periodicity (see Figure 1); IRAS 18566 
-1-0408 is the only (quasi) periodic H2CO maser system known. The flares have a periodicity 
of ~237 days, a duration (at half-maximum) of ~30 days, and an order of magnitude increase 
in flux density. The peak flare intensity and maser intensity during the quiescent state show a 
monotonic decay (Figure 1). Possible causes of the decay include changes in maser pumping 
and gas motions affecting beaming. 

The flares are not strictly periodic. The period derived from the December 2006 and 

November 2008 flares is 238±4 days, and the period derived from the December 2006 and 
March 2008 flares is 234±1 days. The flare of August 2009 was well-sampled during the 
rising part of the flare, but we do not have a good determination of its peak. Assuming 
that the measurement of August 15, 2009 corresponds to the maximum of the flare, then 
the period would be 244±4 days. If the 2002 flare were part of a periodic cycle of outbursts, 
then the period would be 248±5 days (the uncertainties are 3a errors from autocorrelation 
analysis). Most of the recent flares have a shorter period than the one estimated from the 
data prior to 2007. 

The combined effect of (quasi) periodic flares, monotonic flux density decay, and ten- 
dency toward decreasing period as a function of time are exemplifled by the dotted-line in 
Figure 1. The dotted-line was generated as follows: 1) the light-curve from May 2006 to 
January 2007 was used as a template; 2) the template was extrapolated assuming an expo- 
nential decrease in the interval between consecutive flare peaks (r = 237 • 2~^/^^ days, where 
j = 1 for the interval between the December 2006 and August 2007 flares, j = 2 for the 
interval between the August 2007 and March 2008 flares, etc.), and an exponential decay in 
the maser flux density (5*,^ oc 2^*^*^*'^'='^/^^® mJy, where tj-cf is the Julian date of the December 
2006 flare peak). This simple extrapolation (representative fit) based on the 2006 - 2007 
observations roughly reproduces the peak date and variability profile of the fiares, as well 
as the fiux density measurements between 2003 and 2005. The fiare of August 2009 shows 
that the period is neither constant nor strictly monotonically decreasing. Nevertheless, the 
occurrence of the fiares has been regular enough to successfully predict the approximate date 
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of the 2008 and 2009 flares to schedule observations. Indeed, the MERLIN observations were 
scheduled to coincide with the predicted flare of April 2008 (Figure 2). 

We also monitored the 6.7 GHz CH3OH transition with Arecibo. The 6.7 GHz CH3OH 
spectrum shows a group of nine maser components spread between 78 and 88 km s~^. Figure 2 
shows the H2CO and CH3OH spectra obtained on December 14, 2006, and the light curve of 
the H2CO maser and three CH3OH maser components. We discovered that the CH3OH maser 
component at 87.8 km s~^ (hereafter component 9) shows flares that are similar to the H2CO 
maser flares. Specifically, the peaks of all CH3OH (component 9) flares were simultaneous 
(within 10 days) to the peaks of the H2CO flares@ IRAS 18566+0408 is thus the first system 
where (quasi)periodic correlated maser flares of CH3OH and H2CO molecules have been 
detected, and one of only a few CH3OH periodic flare systems known (e.g., Goedhart et al. 
2004, 2009). 

The other CH3OH maser components also show variability but not as well-correlated to 
the H2CO maser (see for example Figure 2, bottom panel). We discuss briefly the implica- 
tions of the variability of the other CH3OH masers in §4, but an in-depth discussion of the 
variability of all CH3OH masers is beyond the scope of this work. 

4. Discussion 

The velocities of the CH3OH maser component 9 and the H2CO maser (Figure 2, insets) 
differ by 8.3kms~^, and the masers are separated by ~2000AU (12 light-days, 0.32") in 
projection (Figure 3). Given that the two masers showed simultaneous flares, the origin 
of the flares is not located in either one of the maser regions, otherwise one component 
would have always flared before the other. The different variability behavior of the different 
CH3OH components shows that the flares are not caused by an homogeneous and large-scale 
change in the background radio continuum, e.g., in the flare of 2006, the CH3OH component 
3 did not show a flare even though it is located in between H2CO and CH3OH masers that 
showed a flare (compare the position of the masers shown in Figure 3 with the light curves of 
Figure 2). We therefore find unlikely that the flares are caused by a change in the background 
radio continuum. 

It is also unlikely that the flares are triggered by a propagating sound/density wave 
or shock front. For example, assuming that the CH3OH component 9 and H2CO masers 
are located at the same distance from the source that triggers the variability (explaining 



^The 2008 March flare is a possible exception. Two CH3OH measurements separated by 14 days gave the 
same flux density within 3ct, which implies an uncertainly on the peak flare of ~14 days. 
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simultaneous flares) , then the minimum distance between the driving source and the masers 
is ~ 1000 AU (i.e., masers and driving source located in the plane of the sky; see Figure 3). 
This minimum distance implies that the travel time of a lOkms^^ wave from its origin to 
the maser regions is > 500 yr, which is far greater than the time scale of the flares (i.e., ~30 
days flare duration, ~237 days periodicity). Faster (J-)shocks (100 km s^^ and greater) are 
unlikely because the molecules would disassociate with the passing of a single shock front 
(e.g., HoUenbach & McKee 1989, Garay et al. 2002), instead of showing periodic variability. 
The absence of significant changes in line peak velocity and linewidth, as well as the smooth 
decay of the maser emission during the quiescent phase argue against the J-shock hypothesis. 

In contrast, a radiative origin of the flares appears more likely. If the two maser regions 
are in the plane of the sky with the source of pumping radiation between them (assumed 
hereafter to be the massive protostar or its surroundings), then the radiation front would 
take only about 6 days to reach the maser regions. This travel time provides a better match 
to the time scales of the flares. In addition, changes in maser gain due to variability of the 
radiation field are consistent with the excitation mechanism of Class II CH3OH masers (i.e., 
infrared pumping, Cragg et al. 2005). 

Given that the CH3OH and H2CO masers show correlated variabihty, the excitation 
mechanism of the two maser species must be similar. The excitation mechanism of H2CO 
masers has been a controversial topic (Araya et al. 2006, 2007c; Hoffman et al. 2003); our 
discovery of correlated H2CO and CH3OH masers shows that the excitation mechanism of 
H2CO masers is probably infrared pumping. 

The periodic behavior of the flares could be caused by a variety of astrophysical pro- 
cesses; chief among these are stellar pulsations and binary induced variability. Some types of 
pulsating variable stars have periods of hundreds of days (e.g., Eyer et al. 2008). Neverthe- 
less, the light curves of such pulsating stars are signiflcantly different from the ffare-like peri- 
odic variability of the H2CO and CH3OH masers. Pulsations of massive protostars (and/or 
inner accretion disks) cannot be ruled out until detailed numerical modeling is conducted. 

An alternative possibility appears more likely: variability modulated by a binary system. 
Binarity and multiplicity are ubiquitous characteristics of massive stars on both theoretical 
and observational grounds (e.g., Krumholz et al. 2009; Kraus et al. 2007). The orbital 
periods of visual and spectroscopic massive binaries range between a few days to many years 
(e.g., Zinnecker & Yorke 2007), encompassing the time scale of the periodic maser flares in 
IRAS 18566-1-0408. In addition, the observed misahgnment between an ionized jet, a Spitzer 
IRAC 4.5 fim excess, and a SiO outflow could be caused by precession due to a binary system 
(Araya et al. 2007c, Zhang et al. 2007). 
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Recently, van der Walt et al. (2009) proposed that the CH3OH flares in G9.62+0.20E 
(which have a periodicity very similar to the H2CO and CH3OH flares in IRAS 18566+0408) 
are modulated by variability from a colliding wind binary, responsible for a change in the 
background radio continuum and/or pumping radiation. Even though such a colliding wind 
model reproduces the variability of the CH3OH maser flares in G9.62+0.20E, we consider the 
model unlikely in the case of IRAS 18566+0408. The H2CO flares are not strictly periodic 
(Figure 1), thus, an additional, stochiastic mechanism is implicated. Minor local turbulence 
could superimpose irregularities on the response to a periodic cause but does not explain the 
combination of decay and irregularity. 

Here we propose an alternative scenario for the maser flares in IRAS 18566+0408: pe- 
riodic accretion of circumbinary disk material. This process has been predicted and has 
observational evidence in the case of some young, low mass binaries (Artymowicz & Lubow 
1996; Jensen et al. 2007; Giinther & Kley 2002; Mundt et al. 2010). In this scenario, mate- 
rial from the circumbinary disk is accreted onto the protostars or accretion disks, heating the 
dust and increasing the infrared radiation held, resulting in higher microwave amplification 
due to greater maser gain. For example, the smoothed particle hydrodynamics simulations 
by Artymowicz & Lubow (1996) show that a binary system with mass ratio 0.79 and ec- 
centricity e = 0.5, will experience bursts of accretion onto the binary components with a 
dimensionless time-dependence that is quite similar to the H2CO light-curve. In a binary 
with mass ratio ~0.8 (e.g., a 20 and 16 M© binary) and e = 0.5, a periodicity of ~240days is 
expected if the semimajor axis of the most massive component is ~1.1 AU. The short flare 
events traced by H2CO and CH3OH masers in IRAS 18566+0408 (this work) and CH3OH 
masers in G9.62+0.20E (Goedhart et al. 2004) could be caused by an orbital conflguration 
similar to the one discussed above. Less flarc-likc (more undulated and/or aperiodic) vari- 
ability seen in other CH3OH maser sources (Goedhart et al. 2004) would be expected from 
binary systems with lower eccentricities (sec Artymowicz & Lubow 1996). This is, the model 
of Artymowicz & Lubow (1996) could explain not only the time scale for the periodicity, but 
also the duration of the flaring events, depending essentially on the eccentricity of the binary 
system. Thus, periodic maser flares could trace the properties of young massive binaries 
that are actively accreting. Moreover, accretion of circumbinary material may cause mass 
equalization in massive binaries. For example, in the model of Artymowicz & Lubow (1996), 
the lower mass component of the binary is the one with a higher accretion rate (see however 
Giinther & Kley 2002). 
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5. Summciry 

We detected the first (quasi)periodic H2CO maser flare system. The maser is coincident 
with a massive protostellar candidate (IRAS 18566+0408). The periodicity of the flares is 
approximately 237 days. We also detected 6.7 GHz CH3OH flares that are correlated to the 
H2CO outbursts. Regardless of whether the H2CO maser flares reported in this work unveil 
tight massive binaries still undergoing accretion or some other process, our discovery shows 
that short time-scale changes (weeks to months) in physical conditions surrounding massive 
protostars are not random, but that underlying (semi)harmonic mechanisms are at work 
during the process of massive star formation. 

This work was partially supported by a Jansky Fellowship of the NRAO. P.H. acknowl- 
edges partial support from NSF grant AST-0908901. We thank the scheduler officers of 
NAIC, NRAO, and MERLIN for their flexibility on scheduling the observations reported 
here, in particular H. Hernandez at Arecibo. We also thank an anonymous referee for com- 
ments and suggestions that signiflcantly improved the manuscript. 
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Observation Date 



Fig. 1. — Light curve of the H2CO maser flux density in IRAS 18566+0408. Data obtained 
with Arecibo, VLA and GBT are shown with different symbols. Measurements before 2006 
are from Araya et al. (2007b). The dotted line is an extrapolation of the 2006-2007 Arecibo 
observations assuming an exponential decay of the flux density and a slow decrease of the 
flare period (see §3). The solid lines are the functions S^^i = 2^'^*^*''°*''/^'^^ mJy, where 
tref = 2454084 JD (the date of the December 2006 flare peak), and Ai is 85 and 9.0 for the 
upper and lower curves, respectively. 
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Fig. 2. — Spectra of H2CO and CH3OH masers obtained with the Arecibo Telescope on 
December 14, 2006 (insets). We show the Arecibo Telescope light curves from May 2006 to 
November 2009 of the H2CO (top panel) and of three CH3OH maser components (middle and 
bottom panels). The H2CO absorption at 85kms~^ originates from the extended molecular 
cloud where the massive (proto)star is located. The dates of the MERLIN observations 
(Figure 3) are marked with vertical lines. 



- 13 - 



Color: [77.5,89.5] km/s 
78 80 82 



o 
o 
o 

CM 



o 

5 



u 

UJ 

□ 



04 12 15.8 

15.7 

15.8 

15.5 

15.4 

15.3 

15.2 

15.1 

15.0 

14.9 
18 59 10. 







84 



86 



88 




H2CO 



-1300 AU 
~7 light-days) 



_L 



_L 



_L 



_L 



_L 



_L 



02 10.01 10.00 09.99 09.98 09.97 09.96 09.95 

RIGHT ASCENSION (J2000) 
Methanol Peak Intensity = 9.3 Jy/b 

Cont. Levels = 12.0 x (-10, 10, 20, 40, 80, 160, 320, 640) mjy/b 



Fig. 3.— Velocity field of 6.7 GHz CH3OH masers in IRAS 18566+0408 obtained with MER- 
LIN {Ogyn = 72 X 42 mas, Position angle = 22°; April 2008). The labelled components cor- 
respond to the spectral features marked in Figure 2 (inset, lower panel), and the variability 
curve of features 3, 8, and 9 are shown in Figure 2 (lower and middle panels). The "©" 
symbol shows the location of the 6 cm II2CO maser {Ogyn = 0.5" x 0.4", P.A. = 11°, VLA 
observations, Araya et al. 2005). The size of the H2CO symbol represents 3a positional 
uncertainty from the fit. 
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Table 1. H2CO Monitoring Observations 
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ZZ. i 


9 9 
Z.Z 


Arecibo 


1 R On+ 
iO ^JCX. 


ZUU4 


9Q Q 

zo.y 


O.D 


VoriJ i 


1 7 Till 


ZUUO 


QR Q 

oD.y 


7 7 
( . ( 


<orlD i 


9P^ A/Tt-nr 

zo iviay 


ZUUD 


io.4 


9 9 
Z.Z 


Arecibo 


ZD iViay 


ZUUD 


1 Q n 
iy.u 


9 7 
Z. / 


AreciDO 


z / iviay 


ZUUD 


1 Q n 
iy.u 


9 /I 
Z.4 


AreciDO 


9^ Tim 

zo jun. 


ZUUD 


iO.O 


9 Q 

z.y 


AreciDO 


90 Till 
ZU JUl. 


ZUUD 


119 
ii.Z 


9 Q 

Z.y 


Arecibo 


Oi Aug. 


ZUUD 


119 
ii.Z 


1 Q 

i .y 


Arecibo 


uo IJCt. 


ZUUD 


Q 

y.o 


1 Q 

i.y 


Arecibo 


98 Oof 
Zo ^JCt. 


ZUUD 


14 1 
i4. i 


9 R 

Z.D 


riieClDO 


09 F>Qn 

uz uec. 


ZUUD 


44.0 


9 9 
Z.Z 


Arecibo 


1^ uec. 


ZUUD 


1 

oO.i 


9 9 
Z.Z 


rireciDO 


ZD uec. 


ZUUD 


/1 1 9 
4i.Z 


9 9 
Z.Z 


AreciDO 


fl/l Ton 

U4 Jan. 


9007 
ZUU / 


91^ 1 
ZO.i 


9 1 
Z.i 


/ireciDO 


11 Jan. 


2007 


16.9 


1.9 


Arecibo 


18 Jan. 


2007 


18.4 


3.1 


Arecibo 


15 Feb. 


2007 


10.7 


2.1 


Arecibo 


04 Mar. 


2007 


9.9 


2.1 


Arecibo 


01 Apr. 


2007 


8.5 


2.5 


Arecibo 


14 Aug. 


2007 


70 


12 


VLA-A 


24 Aug. 


2007 


28.7 


4.8 


VLA-A 


01 Nov. 


2007 


10.4 


1.8 


VLA-B 


13 Jan. 


2008 


4.2 


1.8 


Arecibo 


26 Feb. 


2008 


6.7 


2.2 


Arecibo 
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Table 1 — Continued 



Date (mJy) RMS (mJy) Telescope 



19 Mar. 


2008 


39.8 


2.9 


Arecibo 


22 Mar. 


2008 


53.8 


3.2 


VLA-C 


27 Mar. 


2008 


40.0 


2.5 


Arecibo 


01 Apr. 


2008 


39.6 


3.9 


VLA-C 


01 Apr. 


2008 


34.5 


2.7 


Arecibo 


07 Apr. 


2008 


32.0 


5.3 


GBT 


11 Apr. 


2008 


21.8 


2.2 


Arecibo 


15 Apr. 


2008 


38.0 


7.4 


GBT 


05 May 


2008 


< 15 


5.0 


GBT 


14 May 


2008 


< 15 


5.0 


GBT 


18 May 


2008 


5.1 


2.3 


Arecibo 


23 May 


2008 


7.4 


2.7 


VLA-C 


26 May 


2008 


7.9 


1.5 


VLA-C 


14 Oct. 


2008 


4.8 


1.9 


Arecibo 


22 Oct. 


2008 


< 7.5 


2.5 


GBT 


23 Oct. 


2008 


6.8 


2.0 


Arecibo 


27 Oct. 


2008 


7.9 


2.2 


Arecibo 


02 Nov. 


2008 


8.7 


2.2 


Arecibo 


06 Nov. 


2008 


7.7 


2.4 


Arecibo 


10 Nov. 


2008 


13.3 


2.4 


Arecibo 


14 Nov. 


2008 


28.7 


2.3 


Arecibo 


18 Nov. 


2008 


29.4 


1.7 


Arecibo 


22 Nov. 


2008 


47.3 


2.3 


Arecibo 


02 Dec. 


2008 


27.4 


3.3 


Arecibo 


09 Dec. 


2008 


27.1 


2.6 


Arecibo 


10 Jan. 


2009 


9.4 


3.2 


Arecibo 


11 Mar. 


2009 


5.1 


2.1 


Arecibo 


14 Apr. 


2009 


< 9.3 


3.1 


Arecibo 


15 May 


2009 


< 9.6 


3.2 


Arecibo 


12 Jun. 


2009 


8.8 


2.9 


Arecibo 


27 Jun. 


2009 


7.0 


2.9 


Arecibo 
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Table 1 — Continued 



Date 


Sz/ (mJy) 


RMS (mJy) 


Telescope 


03 Jul. 2009 


11.3 


3.4 


Arccibo 


09 Jul. 2009 


10.1 


2.1 


Arecibo 


17 Jul. 2009 


15.4 


3.7 


Arecibo 


27 Jul. 2009 


20.8 


3.2 


Arecibo 


15 Aug. 2009 


30.0 


2.3 


Arecibo 


30 Sep. 2009 


8.4 


3.0 


Arecibo 


07 Nov. 2009 


6.8 


3.1 


Arecibo 



Note. — Data prior 2006 is from Araya et al. (2004, 
2005, 2007a, 2008). The synthesized beams of our VLA ob- 
servations in A, B, and C configurations were approximately 
0.5", 1.7", and 5", respectively. 



